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ABSTRACT 

We present radio and optical observations of Holmberg I (Ho I) , a member of the M 81 group of galaxies 
(distance ~ 3.6 Mpc). Hoi is a low-mass, low surface-brightness dwarf galaxy. High-resolution multi- 
array Very Large Array observations in the line of neutral hydrogen (H I) reveal a supergiant shell 
(diameter 1.7 kpc) which covers about half the optical extent of Hoi and which comprises 75% of the 
total H I content (total H I mass: 1.1 x 10 s Mq). We estimate the scaleheight of the H I layer to be 
250pc < h < 550pc. We set a tentative upper limit to the dark matter content of < 3.1 x 10 s Mq. 
The H I data are complemented by deep, optical UBV(RI) C observations and narrow band Ha imaging 
obtained at the Calar Alto 2.2m telescope. We find Mhi/Lb — 1-1 Mq/Lbq- The total visible (stars 
plus gas) mass of Hoi adds up to 2.4 x 10 8 A4q. This leads to a total mass of < 5.5 x 10 s A4q and an 
inclination for Hoi of 10° < i < 14°. 

The origin of Hoi's peculiar H I morphology is discussed in terms of a supergiant shell created by 
strong stellar winds and supernova explosions. We estimate that the energy deposited falls in the range 
of 1.2 x 10 53 erg < E < 2.6 x 10 53 erg (equivalent to 120-260 type II SN explosions). From a comparison 
with isochrones as well as from dynamical modeling based on the H I data we derive an age for the 
supergiant H I shell of ~ 80 ± 20Myr. The morphological center of Hoi (i.e., the center of the ring) 
is offset by 0.75 kpc with respect to the dynamical center. Within the interior of the shell the light 
distribution is exponential with a rather shallow gradient and blue optical colors. Beyond a radius 
corresponding to an H I column density of ~ 10 21 cm -2 , the putative star formation threshold, the disk 
becomes considerably redder and the slope for the exponential light distributions steepens. We attribute 
this to a uniform star formation activity in the recent past within the central 2 kpc of Ho I. Color- 
magnitude diagrams based on our CCD data show that the youngest stars, with ages of 15-30 Myr, are 
situated along the inside of the rim of the giant H I shell which is where we also find some faint H II 
regions. It is speculated that these stars are the result of secondary star formation on the rim of the 
shell. Based on the global morphology and velocity dispersion as well as the location of the H II regions 
we find evidence for ram pressure within the M 81 group. Finally, we discuss the likelihood of Ho I having 
lost some of its interstellar material to the intergalactic medium ( "blow-out" scenario) . 

Subject headings: galaxies: individual (Holmberg I, UGC 5139, DDO 63) — galaxies: irregular — 
galaxies: dwarf — galaxies: photometry — ISM: bubbles — ISM: HI 



1. INTRODUCTION 

Star formation (SF), and subsequent evolution of the 
most massive stars, is believed to play a dominant role in 



shaping the interstellar medium (ISM) of galaxies. Stars 
with masses larger than first produce wind-blown 

bubbles within which core-collapse supernovae (SNe) oc- 



1 Visiting Astronomer, German-Spanish Astronomical Centre, Calar Alto, operated by the Max— Planck-Institute for Astronomy, Heidelberg, 
jointly with the Spanish National Commission for Astronomy. 
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cur, typically on time scales of 0.3 - 4 x 10 7 years (Lei- 
therer et al. 1999). Each of these type II SNe releases a 
total energy of <~ (2 — 3) x 10 53 erg, mostly in the form 
of neutrinos of all flavors (for a review see Burrows 2000) . 
Only about 1% (<~ 10 51 erg) of this energy is deposited as 
kinetic energy in the environment. It is this fraction which 
ultimately creates features in the ISM with typical dimen- 
sions of tens of parsecs known as "bubbles" or "shells," 
depending whether one concentrates on the inferred three- 
dimensional morphology or the observed two-dimensional 
shape. As SF tends to occur in associations and clus- 
ters of stars one finds that the combined effects of stellar 
winds and SNe of an ensemble of high-mass stars can cre- 
ate larger features such as "superbubbles" or "supergiant 
shells" with radii of up to more than 1 kpc. Recent mod- 
els of an evolving young association with a total stellar 
mass of 10 4 A^o predict a total kinetic energy input of 
~ 10 53 erg (Leitherer etal. 1999). Once no further energy 
input is provided, i.e., after the lowest mass stars which 
end their lives as supernovae have disappeared, the ex- 
pansion velocity of the SN-driven bubbles and shells will 
decrease to values similar to the velocity dispersion of the 
surrounding ISM. Eventually they will stall. This state of 
affairs is reached sooner for those bubbles which break-out 
of the disk of their host, leaving a morphology of a ring 
rather than a bubble. 

The expanding shock front around each bubble piles 
up the surrounding ISM on its rim (see, e.g., the reviews 
by Tenorio-Tagle & Bodenheimer 1988; Brinks & Walter 
1998, and references therein) and, if densities and temper- 
atures are appropriate, new SF can occur. This process, 
known as "propagating star formation" (e.g., Mueller & 
Arnett 1976; Gerola & Seiden 1978; Elmegreen 1994), ex- 
tends the time span for energy input and can lead to even 
larger shells, provided the bubble has not broken out of the 
disk. The recurrent deposition of such huge amounts of en- 
ergy and their transformation into kinetic energy leads to 
a heating of the ISM, providing a natural energy source 
for the observed velocity dispersion of the (neutral) gas 
component. 

It should be noted, though, that observational evidence 
for a stellar origin of superbubbles or supergiant shells is 
scarce. Only in very few cases are stellar clusters found 
within giant H I holes (Stewart & Walter 2000). In most 
cases these structures do not seem to host remnant stellar 
clusters, as was pointed out by Rhode etal. (1999). 

Observationally, Heiles (1979, 1984) was the first to 
study the H I morphology and dynamics of bubbles and 
shells in the ISM of the Milky Way. This was followed by 
an analogous study of the ISM in M 31 by Brinks & Bajaja 
(1986). They found structures up to 700 pc in diameter. 
Considerably larger shells have been found in dwarf irreg- 
ular galaxies (dlrr). This can be understood as follows: 
the low gravitational potential binds the ISM less tightly 
to the disk. As mentioned by Walter & Brinks (1999), 
the H I velocity dispersion is approximately the same in 
dwarf galaxies as in larger, spiral galaxies (<~ 6 — 9kms~ 1 ). 
Hence the scaleheight in dlrrs is larger than in spirals. This 
holds true in absolute as well as in relative terms. As a 
result, the volume density of the ISM is lower compared to 
grand-design spiral galaxies, resulting in a lower resistance 
against the kinetic energy input of SF regions. Secondly, 



the shells can evolve over a longer period. This is due 
to the larger scaleheight which delays and possibly even 
prevents break-out. Moreover, since dlrrs show mostly 
solid body rotation, there is little shear that may destroy 
shells once they have formed. Examples of objects which 
show large H I shell radii are Holmberg II (Puchc et al. 
1992), IC2574 (Walter & Brinks 1999), IC 10 (Wilcots & 
Miller 1998), DDO 47 (Walter & Brinks 2001), NGC 6822 
(de Blok & Walter 2000), and the Magellanic Clouds 
(Staveley-Smith etal. 1997; Kim etal. 1999). 

To first order at least, a region of massive SF is oblivious 
to its environment, dumping the same amount of energy 
into the ISM whether it finds itself in a large spiral galaxy 
or in a dwarf. Therefore, for progressively lower overall 
galaxy mass, the impact of SF and its aftermath becomes 
proportionally more and more important. Theoretical 
models, e.g., by De Young & Heckman (1994), Mac Low & 
Ferrara (1999), and Ferrara & Tolstoy (2000), predict the 
loss of mctal-cnrichcd interstellar material associated with 
low-mass dwarf galaxies {M. v i s Me ^5 10 9 Mq) in the form 
of shells breaking out of the disk ("blow-out" scenario). 
As dlrrs are dominant in number in the universe and are 
considered to be the building blocks of larger galaxies in 
"bottom-up" scenarios of galaxy formation, this may have 
implications with respect to the enrichment of the inter- 
galactic medium (IGM) at larger look-back times. 

In some low-mass dwarf galaxies the ISM shows an al- 
most perfect ring feature in H I, more or less centered 
on the optical galaxy. According to the models discussed 
above this may be indicative of a recent starburst in the 
center of the galaxy which created the cavity. In these 
galaxies, the total mass may just have been high enough 
to retain the expanding material. Objects with even lower 
mass {M. visible ^ Nfi M.q) would have suffered "blow- 
away" (Ferrara & Tolstoy 2000) of the entire ISM by a 
similar event. Examples of gaseous rings in dwarf irregu- 
lar galaxies which suggest blow-out are Leo A (Young & 
Lo 1996), SagDIG (Young & Lo 1997), M81 dw A (Westp- 
fahl & Puche 1993), Holmberg I (Hoi, Tully etal. 1978), 
Cas 1 (Huchtmeier, priv. comm.), and Sextans A (Skillman 
etal. 1988; Van Dyk, Puche, & Wong 1998). 

Multiwavelength studies, as presented in this paper, are 
indispensable if we wish to understand these unique sys- 
tems. The particular object under study here is Holmberg 
I (Ho I), a dwarf irregular galaxy in the nearby M 81 group 
of galaxies. This object was detected by Holmberg (1950) 
who optically searched for members of the M81 group. 
Since then, relatively few papers have been dedicated to 
detailed studies of Hoi. Sandage & Tammann (1974) and 
Tikhonov etal. (1992) performed photographic photom- 
etry. Both papers concentrate on obtaining a distance 
measurement using the brightest stars method. Hoessel & 
Daniclson (1984), using early CCD imaging, studied the 
stellar population of Ho I in the Gunn gri system. Radial 
surface brightness profiles have been obtained by Brcmnes, 
Binggeli, & Prugnicl (1998) and Makarova (1999), both of 
which observed Ho I in the context of larger surveys of 
dwarf galaxies. The luminosity and size of H II regions, 
as well as the metallicity of Hoi, have been studied in 
detail by Miller & Hodge (1994, 1996). Westerbork H I 
radio synthesis observations were presented by Tully et al. 
(1978) whereas Puche & Westpfahl (1993) presented a pre- 
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liminary map of the VLA data which are reported in full 
in this paper. 

In this paper we aim to present a consistent picture of 
HolmbergI, dealing with the properties and evolution of 
the stars and gas. In particular, we investigate the pos- 
sible influence of the stellar population on the overall gas 
morphology and the star formation history of this galaxy. 
In Sect. 2 wc present optical and H I observations and de- 
scribe the results in Sect. 3. The wealth of data is discussed 
in Sect. 4 and our conclusions are summarized in Sect. 5. 

2. OBSERVATIONS AND DATA REDUCTION 

In this paper we use radio and optical data to study 
the properties of the dlrr galaxy HolmbergI (Hoi, also 
known as UGC5139 or DD0 63). We combine maps in 
the 21-cm line of atomic neutral hydrogen (H I) with 
UBV(RI) C and Ha observations in the optical regime in 
order to study the gas and the stars, respectively. As Ho I 
is a member of the M81 group we assume it to be at the 
same distance as M81, 3.63±0.34Mpc which corresponds 
to to — M = 27.80mag (Freedman etal. 1994). This is 
not too different from a direct measurement for Ho I itself, 
based on the brightest stars method utilizing photographic 
plates (Sandage & Tammann 1974): m — M = 27.63 mag. 
However, using the same technique, Tikhonov et al. (1992) 
found a distance module of m — M = 29.11 mag, which de- 
viates substantially from the others. 

2.1. Optical observations 

Optical images were obtained at the Calar Alto 2 2.2-m 
telescope. The instrument used was a combination of the 
focal reducer CAFOS (f/4.2) and the CCD Site#ld. This 
CCD has a pixel size of 0'.'53 x 0'.'53, a readout noise of 
5.06 e~, and a gain of 2.3 e~/ADU. The Johnson filters B, 
V, and the Johnson-Cousins filter R c were used with in- 
tegration times of 60min, 45min, and 35min respectively. 
These observations were performed in 1999 January under 
seeing conditions of ~ l'/4. In 2000 January we completed 
the UBV(RI) C imaging with the U (100 min integration 
time) and I c (90 min) filters. The seeing was again <~ l'/4. 
Unfortunately the I c exposure suffers from fringing at the 
w 2.5% level of the global sky background. This results 
in the loss of faint structures and a reduced limiting mag- 
nitude. Additionally, an Ha image (1999 January) with 
60 min integration time was obtained. Due to the band- 
width of the filter (12 nm), this image still contains some 
[N II] emission. All images were bias-subtracted using an 
average value of the overscan region. Flatfield exposures 
were produced by taking the median of different skyflats of 
the same night. All reduction steps were performed using 
the IRAF 3 (Tody 1993) software package. 

All observations were split into three or more exposures, 
allowing a straightforward elimination of cosmic ray pix- 

2 The Calar Alto Observatory is operated by the Max-Planck-Institut fur Astronomie (Heidelberg) jointly with the Spanish Comision Nacional 
de Astronomi'a. 

3 IRAF is distributed by National Optical Astronomy Observatories, which are operated by the Association of Universities for Research in 
Astronomy, Inc., under cooperative agreement with the National Science Foundation. 

4 The Very Large Array (VLA) is a telescope of the National Radio Astronomy Observatory (NRAO) which is operated by Associated 
Universities, Inc., under a cooperative agreement with the National Science Foundation. 

5 The Astronomical Image Processing System (AIPS) was developed by the NRAO. 

6 MIRIAD was developed by the ATNF (Australia Telescope National Facility) which is part of the CSIRO (Commonwealth Scientific and 
Industrial Research Organization). 



els and other image defects by applying the CCD clip- 
ping algorithm CCDCLIP. Calibration was performed with 
<~ 15 Landolt standard stars (Landolt 1992) solving for 
the first-order zero point, the airmass, and color terms. 
Fitting of the calibration coefficients was done using the 
task fitparams. The systematic error in the calibration 
is estimated to be 0.05 mag in B, V, R c , and 0.1 mag in U 
and I c . 

For surface brightness measurements, notoriously bright 
foreground stars were replaced by a Gaussian background 
distribution, equal in mean and scatter to the ambient 
brightness. We applied a median filter of 3'.'2 in size in 
order to achieve smoothness. The data were corrected 
for Galactic extinction EB-v=0-048mag at the position 
of Hoi (Schlegel, Finkbcincr, & Davis 1998) which was 
converted to the Johnson(-Cousins) bandpasses following 
Cardelli, Clayton, & Mathis (1989). In the case of surface 
brightness measurements one is interested in diffuse emis- 
sion and cannot apply the aperture used for single stars. 
Therefore we had to compare the total flux F to t of a point 
source with the flux F aper within the aperture used. To 
determine F tot , we opened the aperture for several stars 
in steps of 0V5 and measured the corresponding flux until 
reaching a constant value in the limit. We assumed this 
limit to be F tot . The difference F tot — F aper leads to cor- 
rections of 0.06 mag in B, V, R c , 0.11 in U and 0.30 in 

Ic- 

Crowded field photometry was performed with 
DAOPHOTII (Stetson, Davis, & Crabtree 1987; Stet- 
son 1990) as implemented in IRAF. We used DAOFIND to 
identify candidates of unresolved sources. Subsequently 
we determined a PSF magnitude by a manual selection 
of ~ 20 sufficiently bright point sources in the field. We 
rejected remaining cosmic ray pixels, slightly resolved 
background galaxies and image artifacts by judiciously 
employing the sharpness and x parameters that were de- 
termined during the PSF fitting. The photometric errors 
of the remaining point sources are displayed in Fig. 1 (be- 
cause of the residual fringing we refrained from attempting 
any stellar crowded field photometry in the I c band). 

EDITOR: PLACE FIGURE 1 HERE. 



2.2. H I observations 

Hoi was observed with the NRAO Very Large Array 4 
(VLA) in B-, C-, and D-configuration (1990 July, 1990 
December, and 1991 March). The parameters of the ob- 
servations are listed in Table 1. For the data reduction we 
used two software packages: calibration was done with 
AIPS 5 and imaging with MIRIAD 6 (Sault, Teuben, & 
Wright 1995). To visualize the data, we made extensive 
use of the ATNF package KARMA (Gooch 1995). 
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EDITOR: PLACE TABLE 1 HERE. 

To start with, unreliable visibilities were flagged by vi- 
sual inspection of the raw uv-d&ta. after which we applied 
flux, complex gain, and bandpass calibrations using the 
tasks SETJY, GETJY, CALIB and BPASS. To extract the 
21cm line emission, a linear interpolation of the contin- 
uum in the itw-plane was applied by taking the average 
of line-free channels on each side of the passband. This 
interpolation was subsequently subtracted from all visibil- 
ities (task UVLIN). As a final step in AIPS, the calibrated 
datasets from all configurations were combined with db- 
CON. 

In the imaging process the Mi>-data were Fourier trans- 
formed using both "natural" and "uniform" weighting 
(task invert within MIRIAD). As both weighting schemes 
have some mutually exclusive advantages, it is desirable 
to find the optimum between low noise and high resolu- 
tion. Therefore the visibilities were also weighted with the 
"robust" -weighting algorithm, as originally developed by 
Briggs (1995). Deconvolution of the resulting cube was 
subsequently performed by applying the Clark Clean- 
algorithm (Clark 1980) (tasks CLEAN and restor). The 
"robust" data cube has a synthesized beam of 8'.'2 x 
7'.'0 in size (approximately 140 pc x 120 pc) with a po- 
sition angle of —73° and an rms noise per channel map 
of 1.4mJybeam _1 (equivalent to a brightness tempera- 
ture of 14.8 K) which corresponds to a column density of 
7xl0 19 cm" 2 . 

In order to decide which emission at a faint level is real 
in the data cube, we produced a "master cube" where all 
noise dominated positions were replaced by "blank" val- 
ues. To do this, we first smoothed the data cube to 20" and 
blanked out all positions where the convolved signal drops 
to below 2.5cr rms. Secondly, the resulting cube was exam- 
ined by eye for emission in each channel map. Only regions 
with signal in three or more consecutive channels were con- 
sidered to contain genuine emission. The same areas which 
were blanked in the "master cube" were blanked as well in 
the "robust" and "natural" -weighted cube. These blanked 
cubes were subsequently used for the calculation of the 
flux and the moment maps, i.e., the velocity-integrated 
H I distribution, the velocity field, and the H I velocity 
dispersion map (task moment in MIRIAD). 

The data cube thus obtained is a combination of a 
CLEANcd image and a residual dirty image. Whereas the 
CLEAN-beam is well defined, the effect of the residual dirty 
beam on the output maps is a priori unknown. This makes 
the determination of the flux level in each channel map un- 
certain and a correction needs to be applied. We used the 
method proposed by Jorsater & van Moorsel (1995) (see 
also Walter & Brinks 1999). The total flux S can be com- 
puted by determining the fluxes of the dirty data cube D, 
the residual cube R, and the clean components convolved 
with the clean beam C via the equation 



We note that we would have overestimated our H I flux by 
~ 70 % if we had not applied this correction. 

3. RESULTS 
7 The "Groningen Image Processing System" 



3.1. Global Optical Properties 

Figure 2 shows a true color image of Ho I based on a 
combination of our UBVR C CCD frames. One can clearly 
discern an underlying population of red objects with su- 
perimposed large collections of blue and presumably young 
stars (especially towards the southern part of the galaxy). 

EDITOR: PLACE FIGURE 2 HERE. 

Before attempting photometry on individual stars or 
star clusters we determined apparent and absolute inte- 
grated magnitudes for Ho I (see Table 2) . Our limiting sur- 
face brightness magnitudes are approximately 26 mag/d" 
(U), 27mag/D" (B), 26mag/D" (V), 26mag/D" (R c ), and 
23mag/D" (I c ). 
Taking 

Lb = 1o -O.4(M b -M B0 ) Lb0 (2) 

and a solar absolute blue magnitude of Mb = 5.50 mag 
(Lang 1992) this leads to L B = 1.0 x 1O 8 L B0 - This value 
is corrected for Galactic foreground extinction correspond- 
ing to E B -v=0.048mag (Schlegel et al. 1998) but not 
for any (unknown) internal extinction within Hoi. How- 
ever, since Miller & Hodge (1996) derived a metallicity 
of 12+log(0/H)=7.7 (~8% solar) and extinction is corre- 
lated with metallicity, we expect the internal extinction to 
be negligible. 

EDITOR: PLACE TABLE 2 HERE. 



3.2. HI morphology and dynamics 

3.2.1. HI distribution 

Fig. 3 shows the "robust" weighted H I channel maps. 
The beam size is indicated in the lower left-hand of each 
panel. It is only a few times larger than the seeing of our 
optical observations. The region to the north-west of the 
center appears in nearly all channel maps, suggesting a 
large velocity dispersion. The area near the central posi- 
tion shows virtually no emission throughout the dataset. 
The total H I diameter of Hoi is <~ 5.8 kpc. 

EDITOR: PLACE FIGURE 3 HERE. 

We used the most sensitive "natural" -weighted cube for 
the determination of the spectrum and the flux. The H I 
spectrum of Ho I is shown in Fig. 4. Its shape is remarkably 
well described by a Gaussian with a FWHM of 27.1 kms -1 
and a central velocity of ~ 140 km s -1 . Although it is not 
clear, a priori, why this should be so, this was noticed 
earlier by Lo, Sargent, & Young (1993) who, at some- 
what lower velocity resolution, found similar H I spectra 
in nine intrinsically faint dwarf galaxies. Stil (1999) also 
mentioned that dwarf galaxies with a central H I hole show 
a single peak in their spectra. We tried to model the H I 
distribution of Hoi with the task GALMOD, which is part 
of the software package GIPSY 7 . We started with the ob- 
served properties of Hoi. Changing parameters like the 
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rotation curve and/or the H I column density distribution 
while keeping the inclination and the dispersion fixed, led 
to a bimodal H I spectrum. On the other hand, vary- 
ing the dispersion and/or the inclination, and leaving the 
column density distribution and the rotation curve un- 
touched, kept the Gaussian shape of the spectrum. 

EDITOR: PLACE FIGURE 4 HERE. 

Integrating the flux over all channels yields 36.0 ± 
4.0 Jykms -1 , corresponding to a total H I mass of 1.1 x 
1O 8 A10. This result places Hoi in the intermediate to 
low-mass range of that found for dwarf irregular galaxies. 
Our values compare well with earlier single-dish measure- 
ments by Dickel & Rood (1978) who found 39.7 Jykms" 1 . 
This indicates that we hardly miss any extended emission 
in our interferometric observations. Note that the value 
quoted by Allen & Shostak (1979) of 49.0 Jykms" 1 is off 
by a fairly large margin. We do not have an explanation for 
this difference other than possible confusion with ambient 
neutral (Galactic) gas within their large beam. The ob- 
served H I mass leads to an Aim/ Lb ratio of 1.1 A4q/Lb q ■ 

A prominent depletion of atomic gas is seen in the cen- 
ter of the integrated H I map (Fig. 5), while the bulk of the 
emission is found in a ring. The highest column densities, 
2.0 x 10 21 cm~ 2 , are encountered towards the south-east, 
where the column density is marginally higher than the 
average value in the ring. 

EDITOR: PLACE FIGURE 5 HERE. 

In Fig. 6 we plot the azimuthally-averaged column den- 
sity versus radius, taking the center of the ring as the 
origin (i.e., the position of lowest H I column density as 
indicated by a cross in Fig. 5). As we will explain below, 
this does not coincide with the dynamical center of the 
galaxy. Figure 6 shows that the ring peaks at a radius of 
52", corresponding to a diameter of <~ 1.7kpc, followed by 
a fairly steep decline. From a radius of 125" onwards the 
H I profile levels off. A determination of the H I mass lo- 
cated in the ring by integrating the radial profile down to 
a radius of 125" yields ~ 8 x 10 7 Mq, which is nearly 75% 
of the total H I mass in the system (see also Sect. 4.7). 
The contrast in column density between the central H I 
depression at a 20 oo = 9 h 40 m 30 s , <5 200 o = 71°11'1"8 and 
the ring-like structure is about a factor of 20, with a cen- 
tral, average value of 6 x 10 19 cm~ 2 . In the following, we 
refer to this location of lowest column density within the 
ring as the morphological center and use it as the center 
of the ring itself. The error of this position is only given 
by the beam of the H I observation. However, the ring 
seems to be slightly elongated in the north-south direc- 
tion. Judging from the rim alone, the genuine center of 
the ring-like structure might differ from the point of low- 
est column density by ~ 10". 

EDITOR: PLACE FIGURE 6 HERE. 

3.2.2. H I velocity field 

Fig. 7 shows the intensity-weighted velocity map of Ho I, 
overlaid as contours on the integrated H I map. The 



projected velocity contours shown span a range between 
130 kms" 1 and 150 kms -1 , implying that Hoi is rotating. 
However, the irregularity of the isovelocity contours shows 
that turbulent motions dominate. This is confirmed by 
inspection of position-velocity (pV) cuts along the kine- 
matical major and minor axis and through the center of 
the hole (Fig. 8). The ring, described in the previous sec- 
tion, is not seen to be expanding. However, this might be 
due to the rather face-on view of Ho I. 

EDITOR: PLACE FIGURE 7 HERE. 



EDITOR: PLACE FIGURE 8 HERE. 

We performed a tilted-ring analysis using the H I data 
to derive the dynamics of Ho I. For this purpose, we 
smoothed the H I datacube to 20" resolution and used the 
task ROTCUR (Begeman 1989) in GIPSY. Model velocity 
fields were subsequently developed with the task velfi and 
subtracted from the observed velocity field. This approach 
helps to set restrictions to the parameters determined by 
ROTCUR. The kinematic center as well as the systemic 
velocity for Ho I were found to be well constrained. The 
dynamical center of Hoi is located some 45" (~ 0.75 kpc) 
to the north of the morphological center, at a systemic 
velocity of (141.5 ± l)kms _1 . The dynamical center is 
indicated by the intersection of the major and minor axis 
in Fig. 5 at a 20 oo = 09 h 40 m 31?6, S 2a0 a = 71° 11/45". The 
uncertainty in determining the center is of order 5". 

We were unfortunately not able to get stable results for 
the inclination (as is usually the case for nearly face-on 
galaxies). Not even the construction of model galaxies 
with the task GALMOD covering the inclination-scalchcight 
parameter space helped to constrain its value. The posi- 
tion angle was determined, within the errors, to have a 
constant value of 50°, measured from north through east. 
Rather than to plot the intrinsic rotational velocities we 
present in Fig. 9 a plot of V x sin(i) as a function of radius. 
The rotation curve reaches a maximum near the peak of 
the H I ring beyond which it starts to decline. The rota- 
tion curve flattens out again for R > 1.5 kpc, just where 
the H I ring ends. It should be borne in mind, though, 
that this might be an artifact due to the fact that both 
the position angle and inclination were kept fixed. Our 
results are in agreement with the earlier, low resolution 
data published by Tully etal. (1978). 

EDITOR: PLACE FIGURE 9 HERE. 

We used Fig. 9 to set an upper limit to the inclination 
of Hoi to i < 14°, by simply equating the dynamical mass 
to the sum of the gaseous and stellar mass, the latter one 
based on the optical blue luminosity (Sect. 3.1) and assum- 
ing an M s tars/ Lb = 1-0. This is, of course, a lower limit 
for the total mass (not accounting for molecular gas and 
dark matter; see also Sect. 4.6). 

The H I velocity dispersion shown in Fig. 10 indicates 
that there are two distinct regions in Hoi. In the southern 
part we find values of <~ 9kms~ 1 , similar to the disper- 
sion encountered in other quiet, gas-rich (dwarf) galaxies 
(see e.g., Stil 1999). In the north-western part, however, 
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significantly higher values of ~ 12kms _1 are found. A 
visual inspection of various spectra obtained at different 
positions reveals that the high velocity dispersion is an 
intrinsic property of the H I in Hoi, and not due to mul- 
tiple velocity components blending along the line of sight. 
Table 3 summarizes the H I and optical characteristics of 
Hoi. 

EDITOR: PLACE FIGURE 10 HERE. 



EDITOR: PLACE TABLE 3 HERE. 



4. DISCUSSION 

4.1. The Light Distribution 

In Fig. 1 1 we present the Johnson-Cousins R c band and 
the Ha image, the latter one with superposed contours of 
the integrated H I emission. The extent of the H I is com- 
parable to that of the underlying red low-surface bright- 
ness stellar component. The optical surface brightness dis- 
tribution of the brighter component is asymmetric, being 
higher towards the southern half, more or less coinciding 
with the H I depletion. It is the same region which shows 
the bluest color and hence suggests a relatively young pop- 
ulation (cf. Fig. 2). At a level of ^ B = 25mag/D" (D 25 ) 
this population extends over 220" x 60", or 3.7 kpc x 1 kpc, 
with a position angle of <~ 112°. 

EDITOR: PLACE FIGURE 11 HERE. 

Fig. 12 shows the azimuthally-averaged UBV(RI) C sur- 
face brightness profiles of Hoi centered on the morpho- 
logical center of the H I hole. As an interesting feature 
we note that the U, B, V, and R c band profiles show a 
change in slope at about the D25 diameter, with a steep 
exponential profile at the outer parts of Ho I, and a shal- 
lower exponential disk towards the center. The difference 
in slope between the inner and outer parts becomes less 
pronounced with increasing wavelength. An elliptical av- 
eraging, following the optical isophotes, leads to a similar 
distribution with the kinks at the same radii. 

If it were not for the absence of a bulge component, the 
surface brightness profile would be classified as a Type II 
profile (Freeman 1970). The kink is located at the outer 
edge of the H I ring, where the H I column density falls 
below the empirical SF threshold of 10 21 cm~ 2 (Skillman 
1987), as indicated by the right-most dotted vertical line 
in Fig. 12. We disregard absorption as a possible mech- 
anism for the change in slope. Given the low (8% solar) 
metallicity, dust is expected to play a modest role, if any. 
Besides, the effect becomes more pronounced as one moves 
outward which would require an increasing dust opacity 
with increasing radius, which is unlikely. 

A more acceptable explanation is that SF has been oc- 
curring fairly uniformly within the inner 74" radius. Be- 
yond this radius, current SF is suppressed (column densi- 
ties are below the canonical star formation threshold) . As 
SF has occurred recently, supposedly superimposed on an 
older underlying component, the light within the SF ra- 
dius is still predominantly blue, hence the rather flat and 



blue color distribution. A similar trend has been observed 
in M 81 dwarf A and might therefore be a general feature 
of low-mass dlrr galaxies (Ott et al., in preparation). 

Fig. 13 compares the azimuthally-averaged B— R c color 
with that of the Ha emission. Since the H II regions are 
distributed in a few patches, one should not take the av- 
eraging too literally. But it does seem significant that 
the Ha emission is found towards the outer edge of the 
H I ring. The current star formation rate (SFR) is low: 
Miller & Hodge (1994) derive a total Ha luminosity of 
4.27 x 10 38 erg s _1 for Ho I, which corresponds to a SFR 
of only 0.004A4 Q yr- 1 . 

The fact that star formation is almost absent can be un- 
derstood in the framework of a universal minimum surface 
density threshold in irregular galaxies (Gallagher & Hunter 
1984; Skillman 1987; Taylor et al. 1994; van Zee, Skillman, 
& Salzer 1998; Hunter, Elmegreen, & van Woerden 2001). 
Its value hovers typically around an H I column density 
of 10 21 cm" 2 which corresponds to a mass surface density 
of 8Mq pc~ 2 . There is a debate going on about its exact 
value and its dependence on Hubble type and metallicity 
(e.g., Franco & Cox 1986). For larger systems, Kenni- 
cutt (1998) finds that a sample of 61 normal disk galax- 
ies show Ha emission above a surface density of typically 
£hi+h« = 10 0,5 A^opc -2 . In other words, the threshold 
for disk galaxies corresponds to ~ 4 x 10 20 cm~ 2 . Inspect- 
ing Fig. 11, SF in Ho I occurs where the H I column density 
exceeds ~ 9 x 10 20 cm~ 2 or 7 A4q pc~ 2 . Hence, Hoi man- 
ages to reach values close to the empirical threshold value 
for dwarf galaxies in a few areas only, i.e., coinciding with 
the densest regions in the ring (note that in the case of 
dwarf galaxies one usually ignores the contributions by He 
or H2 to the mass surface density). 

An alternative explanation for a star formation thresh- 
old is that of gravitational instability. In gas-rich, spi- 
ral galaxies a thin gas disk will undergo gravitational col- 
lapse if its gas surface density exceeds a threshold surface 
density. This critical surface density is: T, cr i t = Q 3 3g G 
(Toomre 1964; Kennicutt 1989) where a is a constant of 
order unity, a the 1-dimensional gas dispersion, and G 
is the gravitational constant; n, the epicyclic frequency, 
is defined as: k = 1.41^(1 + f ) 1/2 . If we consider 
only the H I gas parameters (cf. Tab. 3) and use the 
radii R and velocities V given by the rotation curve and 
take into account the derived inclination range at the lo- 
cation of the supergiant H I ring (cf. Sect. 3.2.2), we find 
20A^ Q pc -2 < E crit < 70.M Q pc~ 2 . The observed mass 
surface density never reaches more than 7 MqPC~ 2 , i.e., 
well below what would be needed for SF to take place. 

There are good reasons to argue that a thin disk approx- 
imation is not valid in the case of thick dwarf galaxies and 
alternative approaches are required (Hunter & Plummer 
1996; Hunter, Elmegreen, & Baker 1998; Wilcots & Miller 
1998; Hunter, Elmegreen, & van Woerden 2001). However, 
observed values for S gos , even when corrected for He, lie 
comfortably below the threshold indicating that the gas 
disk is stable against gravitational instability and as a re- 
sult does not experience large scale star formation. Both 
formalisms, a universal threshold and a gravitational in- 
stability criterion, seem to be able to explain the lack of 
current star formation in Hoi; based on the current data 
we find that it is not possible to decide which of the mech- 
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anisms is the dominant one. 

EDITOR: PLACE FIGURE 12 HERE. 



EDITOR: PLACE FIGURE 13 HERE. 

Most of the Ha emission in Hoi lies on the outer rim 
of the large superbubblc. This is in contrast to what has 
been found to be the case for the shells of, e.g., Holm- 
berg II (Puche etal. 1992) and IC 2574 (Walter & Brinks 
1999), where some of the brightest H II regions tend to co- 
incide with the inner rims of the shells. We'd like to note 
that there is an obvious difference, however. In Holm- 
berg II and IC 2574 most of the rings showing Ha emission 
are still expanding whereas the shell in Hoi has virtually 
stalled. The faint Ha seen just to the outside of the rim 
of Hoi's huge H I shell might thus be due to a sufficient 
amount of material being piled up at the front of the shell. 

4.2. Color-Magnitude Diagrams 

The seeing in our CCD images is <~ l'/4 which corre- 
sponds to a linear resolution of 24 pc at the distance of 
Hoi. This means that only the largest open clusters are 
likely to be resolved. As a consequence, we are not able to 
distinguish between most of the small stellar clusters and 
individual stars. 

To study exclusively objects belonging to Hoi it is 
necessary to discriminate them from unrelated fore- 
ground/background sources (Galactic stars/unresolved 
galaxies). This was done assuming that all sources out- 
side a radius of ~ 1.'9, centered on the H I depression 
(cf. Sect. 4.1), arc unrelated objects. The vast majority 
of these objects have colors (U-B)>0.5, (B-V)>0.8, and 
(V— R)>0.6. We assume that objects with the same col- 
ors that fall within the aforementioned radius of <~ 1.'9 
are background or foreground sources, too, and we ex- 
clude them from further analysis. The color-magnitude 
diagrams of the remaining point sources which likely be- 
long to Hoi are shown in Fig. 14 and 15. 

First, we start with assuming that all observed objects 
within Hoi arc unresolved stellar clusters and compare 
their properties with the well-studied cluster population 
in the Large Magellanic Cloud (LMC), using the compila- 
tion of LMC clusters by Bica et al. (1996). In Fig. 14 we 
overlay their sample on the (U— B, B— V) color-color dia- 
gram and the (V, B— V) color-magnitude diagram (CMD) 
of our objects. To make a proper comparison, we adjusted 
the reddening of the LMC cluster compilation to the red- 
dening towards Hoi and scaled the brightness according 
to the distance of Hoi. Adding up the point sources above 
a V magnitude of 22.2 + 1.5 (B— V) — indicated by a solid 
line in Fig. 14 — yields a number of 40 for Ho I and 593 
for the LMC. In other words, about 15 times more clusters 
would have been observed if the LMC were placed at the 
distance of Hoi. Even correcting for the fact that the LMC 
is a factor of 3 more gas rich than Hoi (LMC H I mass: 
3.1 x 10 8 M Q , Luks & Rohlfs 1992), and scaling the num- 
ber of stellar clusters with this factor, leads still to 5 times 
more clusters in the LMC, compared with the number of 
bright point sources in Hoi. The CMD of Hoi most likely 
consists of both stellar clusters and single stars. Therefore 



the cluster abundance of Ho I might be even lower. This 
is supported by the color-color diagram in Fig. 14, where 
in general the LMC clusters with (B— V)~ 0.6 are fainter 
in the U band compared to the point sources of Ho I. 

A reason for the difference might be the ongoing interac- 
tion of the LMC with the Small Magellanic Cloud and the 
Milky Way. Kennicutt etal. (1987) and Bushouse (1987), 
for example, find an enhanced SFR in galaxies with a close 
companion. Hoi's neighbors, the M81 triplet galaxies (cf. 
Sect. 4.5), are at a much larger projected distance for which 
reason the effect of tidal interactions is thought to play a 
less prominent role. Based on our data we can only say 
something about the bright and presumably rather young 
stellar clusters. The number of old clusters in the LMC 
and in Hoi might be the same, after adjusting for the sizes 
of the two galaxies. 

EDITOR: PLACE FIGURE 14 HERE. 

Second, we discuss the consequences of the assumption 
that all point sources are stars. We overlaid the isochrones 
by Bertelli et al. (1994) on the CMDs in Fig. 15. The cho- 
sen metallicity of Z=0.001 is close to the spectroscopically 
derived 12+log(0/H)=7.7 abundance of the brightest H II 
region (Miller & Hodge 1996). Note, that one cannot di- 
rectly compare the metallicity given by the oxygen abun- 
dance with the heavy element abundance described by 
Z. For the Magellanic Clouds, however, Gilmore & Wyse 
(1991) find a nearly solar [0/Fe]~ —0.3. The isochrones 
have been reddened according to Galactic foreground ex- 
tinction (Schlegel et al. 1998) in order to match the uncor- 
rected distribution of stars in the CMDs. There is some 
evidence of the main sequence (MS) in all CMDs, although 
it is not the most populous region in the diagrams. Apari- 
cio & Gallart (1995) found that crowding effects generally 
shift blue stars to the red and red stars to the blue. This 
finding suggests that the MS might be more populated 
than it appears to be in the CMDs. Together with the ex- 
istence of H II regions, this is an indicator for young, not 
very evolved stars. However, the bulk of stars is at least 
older than 30Myr (log t = 7.5). In addition, we do see 
stars in the helium burning (HcB) blue loop phase. Due 
to the described crowding effects, the number of HcB stars 
might be overestimated. 

EDITOR: PLACE FIGURE 15 HERE. 

The filled circles in the (V, U— B) diagram indicate all 
stars within a radius of Of 9, well inside the main H I ring 
feature. These stars appear to follow isochrones between 
15 and 30Myr (log t = 7.2-7.5), which is about half of the 
age derived from the H I properties of the ring (see also 
Sect. 4.4). All these objects are located at a mean radius 
of ~ 660 pc from the morphological center, following the 
ring-shaped H I morphology (cf. boxes in Fig. 16). The 
small time window supports the interpretation for the sur- 
face brightness profiles given in Sect. 4.1. If the supergiant 
H I shell has indeed been created by SF, part of the trig- 
gered population is most likely this population of stars. 

In Fig. 16 the location of the youngest stars is shown. 
These stars were selected to have (U-B)< -0.7, (B-V)< 
—0.1, and V< 23. They neither lie at the H I depression 
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nor at the peak of the H I ring, but rather they are associ- 
ated with the currently visible H I ring. This we interpret 
as more evidence for the H I ring being the likely remnant 
of a superbubble which by its expansion has been the trig- 
ger for subsequent SF and therefore has defined the recent 
SF history of Hoi. 

EDITOR: PLACE FIGURE 16 HERE. 

In contrast, older and redder stars do not show any pre- 
ferred location in Hoi. Van Dyk, Puche, & Wong (1998) 
found a similar behavior in the nearby dlrr Sextans A, 
which shows a comparable H I morphology. In Sextans A 
the HeB blue-loop stars are found mainly inside the H I 
ring-structure. Comparing the stellar photometry of Ho I 
with Sextans A is difficult, because Hoi is at about twice 
the distance. The resulting severe crowding effects can 
only be alleviated, e.g., by the use of space-based tele- 
scopes or adaptive optics. 

4.3. The Thickness of Hoi 

In order to set limits on the thickness of the H I layer of 
Ho I we proceed as follows. We assume that the H I distri- 
bution perpendicular to the disk is Gaussian-shaped; for 
a face-on galaxy this yields 



A, 



HI 



2-7T h nm, 



(3) 



where Nm is the H I column density, h the ler scalehcight 
and nm the H I particle volume density at the midplane of 
the disk. The equation connecting the scaleheight to the 
velocity dispersion is given by van der Kruit (1981): 



h{R) = 



^/4irG Ptot (0,R) 



(4) 



Here, a gas is the one-dimensional H I velocity dispersion, 
G the universal gravitational constant, and p tot (0,i?) the 
total mass density in the disk at radius R and z = 0. Using 
both equations we derive 
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with m p the proton mass. For an average H I column 
density of 3.9 x 10 20 cm -2 , an average velocity disper- 
sion of 9kms -1 , and adopting for ptot/pm a constant 
ratio of 2.2 (assuming p to t/pm = Mtot/Mm] M to t = 
Mm + Mn e + Al stars, with M st ars/L B = 1, a correction 
for the contribution of helium of Mac = 0.3A4m and ig- 
noring any dark matter contribution) , we derive an upper 
limit of ft < 550 pc, corresponding to a FWHM thickness 
of< 1300pc. 

This is comparable to what has been found by Puche et al. 
(1992), Walter & Brinks (1999), and Walter & Brinks 
(2001) for the dlrrs Holmberg II, IC 2574, and DD0 47. 
From this and Eq. 3 we derive an average H I volume den- 
sity of nm ^ 0.10 cm -3 at the midplane of the disk. 



4.4. Evolution of the supergiant H I shell 

AAA. Energy needed to create the shell 

As mentioned in Sect. 1, the standard scenario to ex- 
plain shell structures is the deposition of vast amounts of 
energy via SN explosions of the most massive stars within 
a region of star formation, such as an OB association or 
a stellar cluster. Hoi shows evidence for this scenario as 
we see a concentration of bright, blue stars within the H I 
ring (see Sect. 3.1). The effect of the energy released by 
one SN in the post-Sedov phase on the ambient medium 
has been evaluated by Chevalier (1974). The kinetic en- 
ergy deposited in the ISM can be derived from 



E = 5.3 x 10 



/_no V 

Vcm -3 / 



3.12 



kms 



1.40 



erg, 
(6) 

where no is the volume density of the gas in the ambient 
medium, R the radius, and v the expansion velocity of the 
bubble. 

It has become standard practice to take this equation 
and apply it to the input of an entire stellar cluster. Since 
we do not see the ring expanding, we can assume that 
the shell has more or less stalled at the radius where the 
expansion velocity became comparable to the H I disper- 
sion of the ambient medium (9kms -1 ). However, due to 
the rather low inclination of Ho I, we cannot exclude some 
expansion of the supergiant shell which is not visible in 
the data. As volume density we take the previously calcu- 
lated value for the average H I volume density of 0.10 cm -3 
(Sect. 4.3). Based on Fig. 6 (see also Sect. 3.2.1) we derive 
a radius of the shell of 850 pc. This leads to an energy re- 
quirement of < 1.2 x 10 53 erg to create the supergiant shell, 
which is equivalent to the kinetic energy of ~ 120 type II 
SNe (this value increases by a factor of about two when 
taking into dark matter into account; see Sect. 4.6). 

4.4.2. Age of the shell 

We approximated the evolution of the shell by compar- 
ing it to models of a Sedov expansion phase (Sedov 1959; 
Mac Low & McCray 1988; Ehlerova etal. 1997). The ra- 
dial evolution R(t) of a supergiant shell which is driven by 
a SN rate Asn (each SN provides a kinetic energy E$n) in 
a medium with the particle volume density no and a mean 
molecular mass \i is given by 



R(t) = 53.1 



' N SN \ ( E SN \ 
Myr -1 j V10 51 erg/ 



1/5 
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Myr 
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pc 



(7) 



(note that the expansion velocity goes as V{t) = R{t)). 

If we assume that all of the 120 SNe (Sect. 4.4.1) went 
off within 4 x 10 7 Myr, the longest period for massive stars 
to exist (see, e.g., Leitherer et al. 1999), we derive an aver- 
age SN rate of 3 Myr -1 . Furthermore, we will use again a 
density of 0.10cm -3 (Sect. 4.3), a mean molecular weight 
of 1.3, and an average energy of 10 51 erg per SN. If break- 
out occurs, energy input into the shell becomes ineffec- 
tive, and the structure enters the "snowplow" phase, i.e., 
an expansion with conservation of momentum. McCray & 
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Kafatos (1987) find that for this stage of the evolution of 
a supergiant shell the radius will increase as 

R(t) = R c (8) 

beyond a critical radius R c and a critical time t c , which is 
the radius and the time of the shell right at the moment 
where the shell breaks out of the disk (see also Sect. 4.7). 
Again, the derivative of this equation yields the velocity of 
this phase. Fig. 17 shows the predicted radius and velocity 
at break-out for several critical variables. The line marked 
"SN driven" represents the case without break-out. The 
line shows the steady growth of the shell as a function of 
time. The labels '33,' '24,' etc., represents the expansion 
velocity of the shell at different ages. The family of lines 
marked "snowplow" represent the growth of the shell (at 
a lower rate) after it has broken out of the disk. Note that 
these calculations have been performed for a single shell 
with a central event. Propagating SF and simultaneous 
SF within the hole as suggested by the CMDs and the 
surface brightness distribution (Sect. 4.2 and 4.1) changes 
the energy and the locations of energy input and therefore 
the age calculations (see also Sect. 4.7). 

EDITOR: PLACE FIGURE 17 HERE. 

The age of the supergiant shell can be estimated using 
various approaches: 

• An upper limit for the age of this structure can be 
obtained by simply dividing the radius of the shell 
(0.85 kpc) by the limit of the expansion velocity, 
9kms~ 1 , which leads to 90Myr. 

• The dotted line in Fig. 17 shows the observed 
radius of the shell. Since we found the H I la 
scaleheight to be < 550 pc (Sect. 4.3) we assume 
that the shell broke out of the disk of Hoi, a 
conclusion which is supported by the overall 
morphology (cf. Sect. 4.7). The horizontal dashed 
line of Fig. 17 indicates the size of the shell at 
blow-out. This event should have taken place at 
an age of <~ 15Myr. The intersection of the current 
size of the shell (dotted line) with the evolutionary 
track of the snowplow phase (thick line) leads to a 
total age estimate of ~ 100 Myr. 

• The stars within the H I shell of Ho I and visible in 
the U band are located at a mean radial distance of 
<~ 660 pc (the open squares in Fig. 16). This radius 
was achieved by the shell when it was ~ 35 Myr 
old (see Fig 17, the thicker drawn evolutionary 
track) . Adding the maximum age of these stars 
(30 Myr, see Sect. 4.2) yields an age of ~ 65 Myr. 

• If we assume that the shock front creating the ring 
started as an infinitcsimally thin shock and assume 
that this shock has broadened over time, we find 
that for an equivalent broadening speed of 9kms~ 1 
it would take some 70 Myr to expand isotropically 
by 0.6 kpc (half of the FWHM of the ring). 



• Comparing the integrated light of the remnant of 
the central cluster with the models provided by 
Leitherer etal. (1999) predicts, based on a surface 
color of nu(0) — Mb(0) = —0.41, and adopting a 
metallicity of Z=0.001 and a Salpeter IMF, an age 
of - 60 Myr. 

Given the uncertainties in each of the derived values, we 
conclude an age of the supergiant shell of 80 ± 20 Myr. 
Table 4 summarizes the properties of the H I ring. 

EDITOR: PLACE TABLE 4 HERE. 

4.5. Ram pressure effects 

An obvious aspect of the global optical and H I mor- 
phology is the apparent lopsidedness. We speculate that 
ram pressure may play a role in causing this asymmetry, 
much as was proposed by Bureau & Carignan (2001) in 
the case of Holmbergll. If this is the case, Hoi must be 
falling towards the strongly interacting M 81 triplet con- 
sisting of M 81, M 82, and NGC 3077 (for a description see 
Yun, Ho, & Lo 1994), which lies to the south-east at a pro- 
jected distance about twice the size of the triplet (distance 
Hoi <-> M82: - 120 kpc; M82 NGC 3077: - 80 kpc). 
Ram pressure effects may have compressed the gas in the 
south-eastern part of Ho I, raising it above the critical level 
for SF. It is interesting to note in that context the main 
Ha emission is detected in the same south-eastern quad- 
rant, i.e., the "leading edge" of Hoi. This is somewhat 
similar to what de Boer etal. (1998) found for the LMC. 
On the opposite side, the north-west, one would expect 
higher turbulence (cf. the Coma cluster; Gunn & Gott 
1972). The lincwidth distribution as shown in Fig. 10 cor- 
roborates this scenario. However, it may also be that the 
lopsided distribution of H I in Ho I is due to a (weak) tidal 
pull by the same M81 triplet. We consider detailed calcu- 
lations of either of the two possibilities beyond the scope 
of this paper. 

4.6. The south-east of Ho I 

Besides the big H I hole, Ho I possesses several smaller 
features which are visible in the channel maps (cf. Fig. 3). 
The bulk of these structures are noise-dominated, so we 
refrain from listing all possible cavities. A prominent 
feature, however, is located in the south-eastern part of 
the ring where we find the highest H I column densities. 
The box marked in Fig. 5 is blown up in Fig. 18, where 
we show the H I, the Johnson B-band and Ha emission 
for comparison. This particular H I shell has a diame- 
ter of 0.5 kpc and is filled by blue stars, lending support 
to the picture that massive SF (and the resulting SN ex- 
plosions) is shaping the surrounding ISM. The H I col- 
umn density drops by a factor of ~ 3 towards the center 
of this ring (N m = 6 x 10 20 c m - 2 ; a 200 o = 9 /l 40 m 35 s , 
<520oo = 71°10'00"). A visual inspection of pV cuts at dif- 
ferent position angles centered on the shell, as well as var- 
ious spectra in the natural weighted cube, reveals that the 
shell is still slightly expanding at ~ 6.5 km s -1 . Following 
Sect. 4.4 we derive for a diameter of 0.5 kpc and a density of 
0.10cm -3 an energy requirement of E — 1.7x 10 51 erg, cor- 
responding to only about two supernovae. An upper limit 
to the age of this shell is set by the life expectancy of the 
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least massive star still to go off as supernova, or 40 Myr . 
If we use the same Sedov expansion phase model as in 
Sect. 4.4.2, we estimate for a SN rate of 0.05 Myr -1 an age 
of 15-20 Myr. In contrast to the location of the H II regions 
on the outer rim of the huge supcrbubblc (cf. Sect. 4.1), 
we find Ha regions at the inside edge of the small shell, 
similar to what has been found for other dlrrs (Puche et al. 
1992; Walter & Brinks 1999). For a velocity of 6.5 kms" 1 , 
the rim of the small H I shell needs some 10-15 Myr to pass 
a given point in space. A comparison with the derived age 
might then explain why the H II regions are located at the 
inner rim of this structure. 

The detection of expansion in the z direction of the 
shell implies that it has not broken out of the disk yet. 
The scaleheight of the H I layer must therefore have a 
lower limit of h > 250 pc, assuming an isotropic expan- 
sion. By applying Eq. 5, we find an upper limit for the 
total mass density of p to t ~ 5/9hi within Ho I. If visible 
matter traces the dark matter content, we can estimate a 
total mass for Hoi of < 5.5 x 10 s A4q, which corresponds 
to a limit for the dark matter of < 3.1 x 10 8 Mq. Given 
the relatively small amount involved, this could perhaps 
be accounted for by molecular gas, dust, or even a higher 
■Mstars/ Lb ratio, in which case no dark matter would be 
required to explain the observed rotation curve. Adopting 
our curve leads to an inclination of i > 10°. Employing 
Eq. 3 yields an upper limit for the particle volume density 
of ^0.20 cm~ 3 . We can now use this value to estimate 
an upper limit to the energy calculation of the supcrshcll 
(Sect. 4.4.1): E < 2.6 x 10 53 crg, which is equal to the en- 
ergy input of about 260 SNc. The determination of the 
age of the shell, based on a comparison with models (cf. 
Sect. 4.4.2) does not change substantially; the energy de- 
pends nearly linearly on the volume density in Eq. 6. If the 
density doubles, the energy follows and those two factors 
cancel out in Eq. 7. 

For the small shell, a particle volume density of 
0.20 cm~ 3 in the disk leads to an upper limit of the energy 
needed for its creation of 3.6 x 10 51 erg, which corresponds 
to about 4 SNe. Again, the estimate for its age does not 
change much for the reasons described above. 

EDITOR: PLACE FIGURE 18 HERE. 



4.7. Blow-out 

Whether or not the central big H I hole has suffered 
"blow-out" (i.e., the shell breaking out of the disk and 
losing metal-enriched material to the intergalactic space) 
highly depends on the total mass of the galaxy. In the 
model of Mac Low & Ferrara (1999) and Ferrara & Tol- 
stoy (2000), objects with visible masses below ~ 10 9 .M© 
suffer blow-out. The gas would then be distributed over 
a bigger volume surrounding the former host. They ar- 
gue that for this to work, the velocity of the expanding 
shells would have to exceed the escape velocity of the host 
galaxy. In the case of Ho I we do not see any H I along 
the line of sight to the central hole down to a column den- 
sity of ~ 6 x 10 19 cm~ 2 . The "line of sight" effect may 
contribute to the difference in the measured flux between 
the projected rim of a three dimensional shell and its cen- 
ter. This effect can, however, only account for a change 



of about a factor of two in column density and is insuffi- 
cient to describe the observed H I gradient from the ring 
towards the center (a factor of ~ 20, Fig. 6). 

The calculations performed in Sect. 4.4.2 and displayed 
in Fig. 17 suggest that the shell broke out of the disk, 
since the observed radius of the ring (850 pc) is consid- 
erably larger than the la scaleheight of the gaseous disk 
(< 550 pc). An alternative explanation for the shell to 
stall is that its interior simply has cooled, reducing the 
overpressure. The typical timescale for such a cooling, 
t coo i, is described by McCray & Kafatos (1987): 

tcool = 4C 1 - 5 (A* E 51 ) - 3 n -°' 7 Myr (9) 

where £ corresponds to the metallicity in solar units, A* to 
the number of SN explosions, E 5 i to the energy per SN in 
units of 10 51 erg, and no to the ambient density of the ISM. 
If we use the parameters evaluated for Ho I in Sect. 4.3 and 
Sect. 4.4.1 (A* = 120, n = 0.1 cm -3 ), assume an energy of 
10 51 erg per SN and the Miller & Hodge (1996) metallicity 
of 8% solar, we estimate the cooling time t coo i = 3.7 Gyr. 
This is far too long to play any role in the case of the 
central shell in Hoi. 

If the radius of a shell reaches the la scaleheight of a 
Gaussian shaped halo, the shock will be accelerated per- 
pendicular to the disk, and metal-enriched material will 
be ultimately lost to intergalactic space (Silich & Tenorio- 
Tagle 2001). In this model, neither the pressure of the 
intergalactic medium nor the gravity of the host galaxy 
is able to retain the swept-up gas. This is exactly the 
"blow-out" scenario which was alluded in Sect. 1. 

For radii smaller than the scaleheight h, Gaussian strat- 
ification leads to models quite similar to the ones with uni- 
form density. At R = h the density is 0.68 times the den- 
sity in the midplane. As the radial evolution depends on 
the density as n" 1 / 5 (cf. Eq. 7), the radii and velocities in 
the z direction before break-out will be only slightly larger 
then the ones derived by the Sedov model in Sect. 4.4.2. 
However, if one assumes an exponentially-shaped z dis- 
tribution of the gaseous component, the acceleration and 
therefore blow-out does not happen until a couple of ex- 
ponential scalelengths h exp (Silich & Tenorio-Tagle 1998, 
2001), which is at about 2.4 times the Gaussian scaleheight 
h. In this case, one would expect the shell to be still in 
the Sedov expansion phase with some piled-up gas in the 
z-direction. However, the latter is not supported by the 
observational facts. 

If blow-out indeed has taken place, it is difficult to esti- 
mate how much material was actually lost. The ring seems 
to dominate the galaxy up to a radius of 125" (« 2.1 kpc; 
see Fig. 6). What may the H I distribution of Hoi have 
looked like just before the creation of the supergiant shell? 
We fit different H I profiles to the outer, apparently undis- 
turbed region of the galaxy and extrapolated them back 
towards the center. This was done for a Gaussian, a linear 
and an exponential H I distribution (see Fig. 19). If we 
integrate these profiles, we find that the H I mass corre- 
sponding to the Gaussian fit and the one based on our ob- 
servations are nearly the same, implying that the H I has 
only been moved around. The linear fit underestimates 
the observed H I flux by ~ 30%. The exponential profile 
predicts that ~50% of the H I was lost during the creation 
of the supergiant shell. Although we have no knowledge 
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of the nature of the H I distribution of Ho I before the cre- 
ation of the hole, according to Taylor etal. (1994), many 
dlrrs show an exponential H I distribution. A further ad- 
vantage in adopting such a distribution is that this profile 
exceeds the SF threshold in the inner parts of Ho I out to 
a radius of ~ 1 kpc, which might explain why the burst of 
SF was able to occur offset from the dynamical center, as 
observed. 

EDITOR: PLACE FIGURE 19 HERE. 



5. SUMMARY 

We present high-resolution VLA H I and deep optical 
UBV(RI) C as well as Ha observations of the dlrr galaxy 
Ho I, a member of the M 81 group of galaxies. We find the 
following: 



1. Maps of neutral hydrogen show an impressive 
supcrgiant shell with a diameter of 1.7 kpc. In 
addition, a wealth of small-scale structure is seen 
in the H I data cube. The shell covers about half 
of the H I extent of Hoi (total H I size: <~ 5.8 kpc) 
and contains a substantial fraction (~75%) of 

its total neutral gas content (total H I mass: 
1.1 x 10 8 Mq). We derive the scaleheight of the 
H I layer to be 250 pc < h < 550 pc. 

2. A dynamical analysis shows that turbulence 
plays an important role in Hoi. The velocity 
dispersion in the north-west (~ 12kms _1 ) is 
significantly higher than the overall value of 

~ Okms" 1 . A tilted ring analysis shows that the 
dynamical center is offset by 0.75 kpc from the H I 
morphological center. With the help of a small 
shell, which has not broken out of the disk yet, 
we estimate an upper limit of the dark matter 
content of < 3.1 x 10 8 Mq. This leads to a total 
mass of < 5.5 x 10 8 A^© and an inclination of 
10° < i < 14°. 

3. Extrapolating backwards in time what the original 
H I distribution might have looked like before the 
ring was formed, we estimate a higher neutral gas 
mass than currently observed. Comparing the 
results with published models requires that the 
shell has broken out of the disk. Comparing the 
results with the H I scaleheight and the overall 

H I morphology leads to further circumstantial 
evidence for the occurrence of blow-out in Ho I. 

4. With /x B (0) ~ 24.5mag/D", Hoi belongs to the 
class of low-surface brightness dwarf galaxies. The 
optical extent in the R c band is comparable to the 
H I size. Most of the stellar light falls within the 
H I ring. Azimuthally-averaged surface brightness 
profiles are shallow from the center outwards to a 
radius of 1.2 kpc (where the H I column density 
falls below ~ 10 21 cm" 2 ) and steepen beyond that 
radius. The luminosity profile can be explained 
by a rather uniform star formation history in the 
center of Hoi. 



5. The blue luminosity of Lb = 1.0 x 1O 8 Lb implies 
an Myii/Lb ratio of 1.1 Mq/Lq q . We estimate 

a total visible (stars plus atomic gas) mass of 
2.4xl0 8 7W Q . 

6. H II regions are found on the outer rim of the 
H I shell. This implies that conditions in the 
swept-up H I shell are suitable for star formation 
to commence. Some of the H II regions on the rim 
are centered on a smaller H I shell. 

7. At the distance of Hoi, bright unresolved optical 
sources can be either stellar clusters or single stars. 
A comparison with LMC clusters shows that the 
cluster population near Ho I is much less abundant 
(a factor of < 1/5) than that of the LMC. From 

a comparison with isochrones in the (V, U— B) 
CMD we find that stars within the H I ring are 
relatively young (~ 15-30 Myr). Not surprisingly, 
the youngest stars are found to be close to the H II 
regions on the rim of the H I shell. 

8. The age of the stars within the H I shell as well as 
considerations on the evolution of the gaseous rim 
lead to an age of the shell of - 80 ± 20 Myr. The 
kinetic energy to create the supergiant shell is some 
1.2 x 10 53 erg < E < 2.6 x 10 53 erg (equivalent to 
120-260 type II SN explosions). 

9. The south-eastern side of Hoi is closest to the 
M 81 triplet and shows a steep gradient in the H I 
distribution. The opposite side (north-west) shows 
systematically higher velocity dispersions. Ha 
emission is predominantly found in the south-east. 
All of these features may be indicative of ram 
pressure within the M81 group. Alternatively, 
tidal forces might be shaping the structure of the 
ISM in Hoi at large galactocentric radii. 

Hoi is a fascinating object with an impressive super- 
giant H I shell. Amongst dlrrs with a similar morphology, 
it is one of the very few nearby, low-mass systems which 
can be studied in great detail. From our observations it 
seems likely that the giant H I shell was created by recent 
star formation within the shell (t<~ 80 ± 20 Myr). We can 
only speculate what Ho I will look like in, say, 100 Myr and 
what it looked like before the creation of the H I shell. It 
may well be that Hoi was a blue compact dwarf galaxy 
(BCD) in the recent past, with a centrally peaked H I dis- 
tribution. The active episode of star formation and sub- 
sequent SN explosions may have blown the supergiant H I 
shell that we witness today. In the future the shell may rec- 
ollapse and, in the process, ignite another burst of massive 
star formation ("episodic star formation"). In any case, 
in-depth studies of low-mass and low-metallicity galax- 
ies, such as Hoi, help us understand how the effects of 
massive star formation change the properties of the inter- 
stellar medium. They also give us some insight in to how 
low-mass galaxies may have contributed to the chemical 
enrichment of the intergalactic medium when the universe 
was young. 
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FIG. 1. — Formal photometric errors as a function of magnitude in the UBV(RI) C bands (output of ALLSTAR in daophot). 
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Fig. 3. — H I channel maps of Hoi. The shape of the beam is plotted in the lower left, and the (heliocentric) velocity in the upper left 
corner in units of kms~ 1 (1' corresponds to ~1 kpc). Note the central H I depression in all channel maps and the fact that there is emission 
in almost every channel map in the north— west. Structures at smaller scales, such as rings or holes, are usually visible in two consecutive 
maps only. 
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Fig. 3. — continued 




Fig. 4. — H I spectrum of Hoi. The shape is a nearly perfect Gaussian, as illustrated by the dashed curve, which is a least-squares fit to 
the data points. Velocities are heliocentric. 
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Fig. 5. — Upper panel: The integrated H I emission of Hoi. The beam of 8"2 X 7'.'0 in size is shown in the lower left corner. Lower panel: 
The same figure as above, but with the kincmatical major and minor axes overlaid. The cross indicates the morphological center, i.e., the 
point where the H I density reaches its minimum. The box shows the region discussed in Sect. 4.6. A blow— up of this area is shown in Fig. 18. 




Radius [ " ] 

Fig. 6. — Azimuthally-averagcd H I distribution, taking the morphological center (cf. Fig. 5) as the origin (1' corresponds to ~lkpc). 
The vertical lines represent the radii where the H I column density reaches values sufficiently high for star formation to commence (-/Vjji = 
10 21 cm" 2 ). 
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Fig. 8. — Position— velocity cuts parallel to the major and minor axis and centered on the morphological center of lowest H I column density 
in the hole (cf. Fig. 5). Positive offsets point to the north-west (major axis) and to the south-west (minor axis). 1' corresponds to ~lkpc. 




50 100 150 

Radial offset from dynamical center | *' | 

Fig. 9. — Observed rotational velocity (V X sini) curve of Hoi. The error bars reflect the difference in rotational velocity between the 
approaching and receding sides of the galaxy. 
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Fig. 11. — Upper panel: Optical R c band image of Hoi. Lower panel: Continuum-subtracted Ha image, at the same scale, with contours 
of the integrated H I column density smoothed to 15". Contour levels are drawn at 2 X 10 20 cm~ 2 intervals, starting at 1 X 10 20 cm -2 . The 
optical radiation is mainly emitted in the southern region of Hoi, within a roughly rectangular area which shows an abundance of bright, 
mainly blue stars. The extent over which the R c band emission and the H I distribution are visible are nearly the same. 




Fig. 12. — Azimuthally— averaged UBV(RI) C surface brightness distributions with respect to the morphological center of the H I column 
density distribution (cf. Fig. 6). At a radius of ~ 70" (~ 1.2 kpc) there is a change in slope, diminishing in importance with increasing 
wavelength. The vertical lines correspond to the location of the maximum surface density of the H I ring (solid line) and the range where the 
H I column density exceeds the putative star formation threshold, TVhi = 10 21 cm -2 (dotted lines); see also Figs. 6 and 13. 
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Fig. 13. — Azimuthally-averaged B— R c color (in mag/d", upper graph) and the smoothed Ha emission as a function of radius in arbitrary 
units (lower graph). Note that the Ho peak and the change in slope in the B— R c graph are coincident with the change in slope of all colors 
at a radius of ~ 70" (~ 1.2 kpc); see also Fig. 12. The solid vertical line shows the peak position of the H I ring, whereas the dotted lines are 
drawn at the star formation threshold (ATjji = 10 21 cm -2 ). 
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Fig. 14. — Filled circles: (U— B, B— V) color-color diagram and (V, B— V) color— magnitude diagram of Hoi. Open circles: LMC cluster 
distribution (Bica et al. 1996). A reddening correction corresponding to that towards the position of Hoi has been applied. The calibration 
error in the color— color diagram is plotted in the lower left. In the same diagram, the reddening vector of Hoi is shown in the upper right. 
The errors in the CMD are the same as in Fig. 15. 
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Fig. 15. — Color— magnitude diagrams of Hoi. The errors on the right side of the panels are the magnitude-dependent formal errors given 
by the task ALLSTAR. In the upper left corners, the photometric calibration errors are shown. Isochrones from Bertelli et al. (1994) for a 
mctallicity of Z=0.001 are overplotted with the logarithmic age in years. These isochrones are reddened according to the Galactic foreground 
values at the position of Hoi. The filled circles in the (V, U— B) diagram arc all the stars within a radius of 0'9, inside the H I ring. 
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Fig. 16. — Integrated Hi of Hoi with Ha emission as contours. The markers show the location of young stars selected on the basis of the 
(V, U— B) diagram under the condition (U— B)< —0.7 and V< 23 (crosses). The circles are stars taken from the (V, B— V) diagram with 
(B— V)< —0.1 and V< 23. The open squares show the location of the stars plotted in Fig. 15 as filled circles. 




Fig. 17. — Evolution of the supcrgiant H I shell which dominates Hoi according to the models discussed in the text (Sect. 4.4). After a 
certain period of SN-drivcn expansion the shell enters the "snowplow" phase, due to break-out or cooling. The turn— off points are plotted for 
critical times t c =5-40Myr in steps of 5Myr. The dotted line indicates the observed present-day radius of the supergiant shell. The dashed 
line represents the lcr scaleheight of the H I layer of Hoi. Expansion velocities, derivatives of the radial evolution, are given at the time of 
the break-out. 
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Fig. 18. — Blow-up of the box marked in Fig. 5. Upper panel: Ha emission overlaid as white contours on the H I content. Lower panel: 
Johnson B image. Both images are at the same scale. The ellipse indicates the size and orientation of the H I shell. 




Fig. 19. — Radial H I distribution of Hoi (solid) with three different extrapolations towards the center, in an attempt to estimate the 
amount of H I which might have been present before the energetic events which created the giant hole; exponential (long dashed), Gaussian 
(short dashed) and linear (dotted) extrapolation. The horizontal line labelled "SFT" displays the empirical SF threshold at JVhi = 10 21 cm -2 . 
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VLA configuration 




BCD 


Date of observation 


1990 Jul 23 1990 Dec 02 1991 Mar 07 


Time on source 


340 min 115min 97min 


Secondary Calibrator 


0917+624 0836+710 0836+710 


Coordinates of pointing center (B1950.0) 


a = 09 h 36 m 00 s , 5 = 71°24'47" 


Primary Calibrator 


3C286 


Number of channels 


102 


Velocity range 


321.39 - SS^kms- 1 


Channel width 


2.6kms" 1 


Primary Beam (HPBW) 


32' 


Synthesized Beam (HPBW) 




"natural" weighting 


ll'.'8xll'.'0 


"uniform" weighting 


6'.'0x4'.'7 


"robust" weighting 


8'.'2x7'.'0 


la rms per channel map ("robust" weighting) 


1.40 mjy beam -1 



Table 1 

Properties of the VLA H I observations. 



Johnson (Cousins) Band 


U 


B 


V 


R c 


Ic 


Apparent magnitude [magj 
Absolute magnitude [mag] 
Mean surface brightness [mag/d"] 
Central surface brightness [mag/D"] 


13.23 
-14.57 
25.35 
24.06 


13.25 
-14.55 
25.37 
24.47 


12.49 
-15.31 
24.61 
24.15 


12.40 
-15.40 
24.52 
23.84 


11.87 
-15.93 
23.99 
22.60 


Blue luminosity [£b ] 
Mm/L B [Mq/L Bq ] 


1.0 x 10 s 
1.1 



Table 2 

General optical properties of Ho I. The mean surface brightness is based on a circular area with a radius of 

150". 
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Object Name 


Hoi (UCC5139, DD0 63) 


Position (J2000) 


a = 09' l 40 m 32 s 3 S = +71°10'56" 


Morphological type a 


IAB(s)m 


Adopted distance 


3.6 Mpc (m-M=27.80mag) 
V A 1 kpc 


Scale at this distance 


H I diameter 


5.8 kpc 


Systemic velocity 


141.5kms _1 


Inclination 


10° <i< 14° 


H I linewidth (FWHM) 


27.1 kms" 1 


Total H I flux 


36.0 ± 4.0 Jykms- 1 


Total H I mass 


LlxlO 8 Al0 


Mean H I column density 


3.9xl0 20 cm" 2 


Peak H I column density 


2.0xl0 21 cm- 2 


Average midplane H I particle volume density 


0.10 cm- 3 < n < 0.20 cm- 3 


Average H I la scaleheight 


250 pc <h< 550 pc 


Mean H I velocity dispersion 


9kms _1 


Total mass 


< 5.5 x 10 8 M Q 


Dark matter content 


< 3.1 x 10 s A4q 


Dynamical center (J2000) 


a = 09 /l 40 m 31.6 s , S = 71°11'45" 


Morphological center (J2000) 


a = 9 h 40 m 30 s , 5 = 71°ll'l'.'8 


Optical diameter at = 25 mag/D" 


3.7 kpc x lkpc 


Adopted extinction 13 


E B -v=0.05mag 


Apparent blue luminosity 


13.25 mag 


Absolute blue luminosity 


-14.55 mag = 1.0 x 10 8 L Bq 


Central blue surface brightness 


24.47 mag/D" 


Mean blue surface brightness 


25.69 mag/D" 


Oxygen abundancc c 12 + log(§) 


7.7 


Current star formation ratc d 


0.004 Meyr- 1 



Table 3 

Summarized general properties of Ho I. 



"from the Nearby Extragalactic Database (NED) 
b Schlegel ct al. (1998) 
c Miller & Hodge (1996) 
d Miller & Hodge (1994) 



Radius 


0.85 kpc 


FWHM width 


1.26 kpc 


H I mass 


- 8 x 1O 7 A^ 


Central H I column density 


6 x 10 19 cm- 2 


Average H I column density on the rim 


1.2 x 10 21 cm- 2 


Peak H I column density on the rim 


2.0 x 10 21 cm- 2 


Age 


80 ± 20 Myr 


Energy 


1.2 x 10 53 erg (= 120 SNe) < E 
< 2.6 x 10 53 erg (= 260 SNe) 



Table 4 

Properties of the central supergiant H I shell. 



